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^ ■ Abstract 



X 

5t , We derive the major characteristics of inclusive and diffr active final states 

in deep-inelastic scattering off heavy nuclei for the-high energy (small-x) kine- 
matics in which the limit of complete absorption is reached for the domi- 
nant hadronic fluctuations in the virtual photon (the black-body limit of the 
process). Both the longitudinal and transverse distributions of the leading 
hadrons are found to be strikingly different from the corresponding ones within 
the leading-twist approximation, and hence provide unambiguous signals for 
the onset of the black-body limit. 
The black-body limit (BBL) of deep-inelastic scattering (DIS) from a heavy nucleus 
was first considered by V. Gribov in 1969 before the discovery of perturbative QCD 



(PQCD). He assumed that at ultrahigh energies each hadronic configuration in the photon 
scatters with the same strength, and is absorbed by the nucleus. If the (virtual) photon is 
expressed as a linear superposition of partonic components (characterised by the transverse 
and longitudinal momenta of the individual partons) then the interaction with the target 
does not mix up these partonic components, i.e., they are eigenstates of the scattering matrix. 
This implies that the interaction is diagonal for each hadronic fiuctuation, and features such 
as transverse sizes, invariant masses, and longitudinal momentum distributions are preserved 
by the BBL interaction. Away from the BBL, the diagonality is grossly violated (see, e.g.. 



The proton structure function, -Ff (x, Q"^), measured in DIS at DESY ep collider (HERA) 
exhibits a rapid rise with the energy of the photon-proton subprocess, W, at high energies, or 
small X = —q^l {2Pq), where P is the proton's momentum, and q, = —Q^ are the photon's 
momentum and virtuality. The standard leading-twist approximation (LTA) of PQCD de- 
scribes F2 extremely successfully in terms of a convolution of perturbatively calculable hard 
coefficient functions with quark and gluon parton distributions, which evolve in scale ac- 
cording to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 0] evolution equations. 
The rise of is translated into a rapid rise of the gluon and sea-quark distributions of the 
proton with 1/x at small x. To leading-log accuracy in Q^, this leads to a numerically large 
density of partons at small x, which calls into question the applicability of the LTA at even 
smaller x (for recent reviews, see [§|-0]). This concern is best understood in the target rest 
frame by considering the inelastic interaction cross sections for given partonic configurations 
of the virtual photon. In particular, the expression for the scattering of a small qq dipole of 
transverse diameter dj_ at small x is given, in terms of the leading-log gluon distribution of 
the nuclear target xg"^, by: 



in which the scale = \/d\, with A a logarithmic function of Q^. For large gluon densities, 
this leads to a confiict with unitarity of the 5* matrix for the interaction of spatially small 
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4qcd(4,^) = -dlas{Q^)xg\x,Q') 



(1) 



wave packets in the virtual photon which are almost completely absorbed ^ (especially 
in the case of heavy nuclei, where the gluon density is enhanced by nucleon number A). 
This may be conveniently expressed in terms of the geometrical limit for the inelastic "small 
dipole - nucleus" interaction of 

corresponding to the BBL for the total cross section of ^ttRI^j.^^^. 

Numerical studies performed using Eqs. (|l] and |]) have demonstrated that interactions at 
the upper range of HERA energies in gluon-induced processes may be fairly close to the BBL 
in a wide range of impact parameters (see, e.g., |]8|,p!0|). This is consistent with an analysis 



ri|] of current fits to HERA diffractive data, which lead to diffractive gluon densities which 
are significantly larger than the quark ones (see also ||12|| ). 

For heavy nuclei and central impact parameters, the perturbative x, region, in which 
the interaction should be close to the BBL in both quark and gluon channels, is much 
broader than for nucleons, where scattering amplitudes at peripheral impact parameters are 
important and are far from the BBL. For example, for DIS on Pb at a; ~ 10^'^, at small 
impact parameters the BBL may be reached for < 7(15)GeV^ for the quark (gluon) 
channel , that is, for the region where InQ^/Ag^^ > lna;o/a; (a;o ~ 0.05 -^ 0.1 is the starting 
scale of the evolution at small x). This indicates that DGLAP becomes inapplicable in an 
X, range, where the expansion of a small dipole in the transverse direction (diffusion) 
due to gluon bremsstrahlung is still small and hence the coupling constant is small. This 
indicates an exciting possibility of the existence of the BBL in QCD in the perturbative 
domain for the kinematics which will be probed at CERN Large Hadron Collider (LHC), 
and which could also be probed at Relativistic-Heavy-Ion-Collider (RHIC) (at BNL) in pA 
collisions and at HERA/THERA in eA mode. 

The experimental and theoretical challenge is to find unambiguous signals for the BBL. 
One problem is that predictions for the inclusive structure functions do not contain "smoking 
gun" signatures since they must account for leading twist shadowing, uncertainties in the 



input gluon densities, etc. This is especially acute in the case of the nucleon scattering 
where, even in the BBL, the structure function -Ff Q"^) is expected to increase rather 
rapidly, oc In^l/x, of which In^l/x is due to a Froissart-like increase of the radius 
of interaction r^^^ oc Inl/x (Ogg oc Inl/x) and the remaining Inl/x reflects the infinite 
renormalization of the electric charge [|l|. As discussed in [|10|, the latter results from the 



fact that the photon couples in a pointlike manner to quarks. This also implies that for any 
energy there will always be small configurations in the photon wave function which have 
not yet reached the BBL. However, we expect that the BBL formulae will give a reasonable 
description of structure functions and hard diffractive processes for heavy nuclear target, 
where the blackness of the interaction is enhanced by the number of nucleons. The growth 
in the radius is essentially screened by the nuclear environment so that, for large nuclei, 
F2 oc Inl/x. 

To deduce predictions for the properties of the final states it is useful to work within 
the qq dipole approximation for the photon wave function. In the BBL the elastic and 
inelastic contributions are equal so that the total cross section is given by twice the inelastic 
contribution (summed over all dipole configurations): 



= 2 / |^^f(.,p^,g^)|Vi?^ . (3) 



Here z is the fraction of the photon plus momentum, = qo + q^, carried by the quark, 
z = p^/q^, and p± is the modulus of its transverse momentum in the plane perpendicular 
to q (the antiquark has (1 — z) and —p±). In Eq. (^ we implicitly use the fact that the 
interaction is diagonal in z and p± and black for each configuration. It is convenient to define 

M' = ipl + ml)/izil-z)) (4) 

for the mass squared of the qq system and angle 6 between the direction of the momentum 
of the quark in the centre of mass frame and the photon direction (the transverse plane is 
defined to be perpendicular to this). Neglecting compared to M^, gives 



sin e = 2p±/\/IP, z = {l + cos e)/2 . (5) 



Hence for the nuclear structure functions we get 



Integration over 6 gives the transverse and longitudinal structure functions in the BBL as an 
integral over produced masses. If one takes into account only the suppression induced by the 
square of the nuclear form factor in the rescattering amplitude, then M^g^^ < W'^ /{miyRA) 
and Eqs. and (J^) coincide with the original result of lU. However, in QCD, due to 
color transparency, M^g^^ is determined by the unitarity constraint and is substantially 
smaller than W'^/{mNRA) (for a recent discussion, see The contribution of larger 
is suppressed by a smaller interaction cross section (oc in the LTA with small enough 

gluon densities). In order to visualize manifestations of the BBL, in the following discussion 
we neglect this contribution, which has a stronger A dependence than the BBL term. The 
inclusion of p{M^) corrects the dipole formula, for which Eq. holds, for higher order 
corrections (in a^) contributing at a given M^. A relation between p(M^) and nuclear 
structure functions was also suggested in [[T^ . 

We also obtain, for the ratio of the deeply virtual Compton scattering and elastic ampli- 
tudes, 

M(7^A^7A) _ Jo (M2+(32) 



(M2+Q2)2 

The ability to neglect nondiagonal transitions in the BBL [|I| justifies the removal of the 
integral over the masses in Eqs. d) and (0) and, since diffraction is 50% of the total cross 
section, we immediately find for the spectrum of diffractive masses: 



d cos edM^ d cos edM^ 



(10) 



where e is the photon polarization. Thus, the spectrum of hadrons in the centre of mass 
of the diffractively produced system should be the same as in e^e~ annihilation. Hence, 
the dominant diffractively produced final state will have two jets (with fractal substructure) 
with a distribution in the centre of mass emission angle proportional to 1 + cos^ 9 for the 
transverse case and sin^ for the longitudinal case. The diffractive cross section, integrated 
over 6', is obtained from Eq. (p^) by removing the |(1 + cos^6'), | sin^ ^ factors. It follows 
from Eq. (p!0|) that in the BBL diffractive production of high jets is oc (while in the 
LTA it is oc InQ^) and hence is enhanced: {p±{jet))j, = 3 {p]_{jet))^ = . 

The relative rate and distribution of jet variables for three jet events (originating from 
qqg configurations) will also be the same as in e^e~ annihilation and hence is given by the 
standard expressions for the process e+e^ — > qqg. In addition, in the BBL, the production 
of jets is flavor democratic (weighted by quark charges but unrelated to the quark content 
of the target). 

An important advantage of the diffractive BBL signal is that these features of the diffrac- 
tive final state should hold for < Qbbl even for > Q^bl ^ ^qcd because configu- 
rations with transverse momenta < (5bbl/2 are perturbative but still interact in the black 
regime (and correspond to transverse size fluctuations for which the interaction is already 
black) . 

Another interesting feature of the BBL is the spectrum of leading hadrons in the photon 
fragmentation region. It is essentially given by Eq. (p!0|). Since the distributions in z (or 
equivalently 6) do not depend on M, the jet distribution in z is given by 



d{aT + tOL) 1 + (2^ - 1 



i2 



Exclusive vector meson production in the BBL corresponds, in a sense, to a resurrection 
of the original vector meson dominance model without off-diagonal transitions. The 



amplitude for the vector meson-nucleus interaction is proportion to 2ttR\ (since each config- 
uration in the virtual photon interacts with the same BBL cross section). This is markedly 
different from the requirements |^ for matching the generalised vector dominance model 
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(see, e.g., 0) with QCD in the scahng hmit, where the off-diagonal matrix elements are 
large and lead to strong cancellations. We can factorize out the universal black interaction 
cross section for the dipole interaction from the overlap integral between wave functions of 
the virtual photon and vector mesons to find, for the dominant electroproduction of vector 
mesons, 

Jt ~ It ~ 167r a(M2+Q2)2- _^^2^ ' v-y 

where Ty is the electronic decay width V e^e~, a is the fine-structure constant. Thus 
the parameter-free prediction is that, in the BBL (complete absorption) at large Q^, vector 
meson production cross sections have a l/Q^ behaviour. This is in stark contrast to an 
asymptotic behaviour of l/Q^ predicted in PQCD ||15[ since a factor l/Q^, due to the square 
of the cross section of interaction of a small dipole with the target (color transparency), 
disappears in the BBL. 

In the LTA, the factorization theorem is valid and leads to a universal (i.e., target- 
independent) spectrum of leading particles for scattering off partons of the same flavor. 
Fundamentally, this can be explained by the fact that, in the Breit frame, the fast parton 
which is hit by the photon carries practically all of its light-cone momentum {z 1). As a 
result of QCD evolution, this parton acquires virtuality, ~ Q^, and a rather large transverse 
momentum, kt (which is still ^ Q"^). So, in PQCD, quarks and gluons emitted in the process 
of QCD evolution and in the fragmentation of heavily virtual partons together still carry all 
the photon momentum. In contrast, in the BBL, the leading particles originate from coherent 
diffraction (peripheral collisions) and central highly inelastic collisions. These contributions 
come from the fragmentation of a highly virtual qq pair with similar light-cone fractions of 
longitudinal momenta and large relative transverse momenta [see, e.g., Eqs. ([T0| ) and (|11])]. 

The inclusive spectrum of leading hadrons can be assumed, neglecting energy losses, as 
being due to the independent fragmentation of quark and antiquark of virtualities > Q^, 
with z and p± distributions given by Eqs. (|TD|) and (|TTp (cf. diffractive production of jets 
discussed above). Note that energy losses of partons calculated in the limit of small nuclear 
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parton densities do not lead to a change of the the z fraction carried by a parton, and hence 
do not violate the LTA (see, e.g., [0]). In the BBL, energy losses may be larger, further 
suppressing the spectrum as compared to Eq. (plsl) . 

The independence of fragmentation is justified because large transverse momenta of 
quarks dominate in the photon wave function [cf. Eqs. (0)-(0)] and because of the weakness 
of the final-state interaction between q and g, since is small and the rapidity interval is of 
the order of 1. Obviously, this leads to a gross depletion of the leading hadron spectrum as 
compared to the LTA situation in which leading hadrons are produced in the fragmentation 
region of the parton which carries essentially all momentum of the virtual photon. If we ne- 
glect gluon emissions in the photon wave function, we find, for instance, for the differential 
multiplicity of leading hadrons, (IN^tI^ j dz^ produced by transverse virtual photons, in the 
BBL, 

dz 

Here D''/^[z/y, Q"^) is the fragmentation function of a quark, with any flavor g, into hadrons. 
To simplify Eq. (^) we used = _D^/^ and neglected a rather small difference between 
D""/^ and D'^"/^. 

Calculations using Eq. (|l^) show a strong suppression of the leading hadron spectrum 
as compared to the LTA predictions (see Fig. 1). Moreover we expect a further softening as 

increases, resulting from increased parton emission in the virtual photon wave function: 
progressively more configurations contain extra hard gluons, each fragmenting independently 
in the BBL, further amplifying deviations from the standard LTA predictions. 

Another important signature of the BBL is the hardening of p± distributions with de- 
creasing X (at fixed Q^). Hence, an efficient experimental strategy would be to select leading 
jets in the current fragmentation region and examine their z and p^-dependence as a func- 
tion of X. Qualitatively, the effect of broadening of p<^ distributions is similar to the increase 
of the p_L distribution in , although final states in DIS were not discussed in this model. 



2 j' D'^'\z/y, Q')\{1 + {2y - lf)dy . (13) 



However, due to the distribution of the available p±_ between the produced hadrons, this 
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effect would be more difficult to observe than the z depletion one, and one would require 
detailed modelling. 

An important advantage of nondiffractive scattering is the ability to select scattering at 
central impact parameters (for example, via studies of inelastic interactions as a function 
of the number of nucleons produced in the nucleus fragmentation region). Such a selection 
allows the effective thickness of the nucleus to be increased, as compared to the inclusive 
situation, by a factor of ~ 1.5, and hence allows the BBL to be reached at significantly larger 
X. An important signature will be a strong depletion of the leading hadron spectrum with 
centrality (in the BBL) and a lack of correlations with centrality (in the LTA). 

To summarize, we have demonstrated that the study of final states in DIS off nuclei 
would provide several stringent signals for the onset of the BBL. In contrast to inclusive 
observables they will not depend on details of the input parton distributions, and hence 
will be far more reliable. In the nucleon target case, despite the fact that many photon 
configurations are far away from the BBL, an examination of these final-states observables 
appears to be a promising way to search for precursors of the BBL. In particular, looking 
for a strong depletion of hadrons in the current fragmentation region, in combination with 
the detection of particles in the proton fragmentation region, looks promising. 

We thank A.H. Mueller for a useful discussion and GIF, ARC, PPARC, and DOE for 
support. LF and MS thank INT (University of Washington) for hospitality during the time 
this work was completed. 
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